Electronic Structure of KFe2Se2 from First-Principles Calculations 
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Electronic structure and magnetic properties for iron-selenide KFe2Se2 are studied by first- 
principles calculations. The ground state is coUinear antiferromagnetic with calculated 2.26 
magnetic moment on Fe atoms; and the Ji, J2 coupling strengths are calculated to be 0.038 eV and 
0.029 eV. The states around Ef are dominated by the Fe-3d orbitals which hybridize noticeably 
to the Se-4p orbitals. While the band structure of KFe2Se2 is similar to a heavily electron-doped 
BaFe2As2 or FeSe system, the Fermi surface of KFe2Se2 is much closer to FeSe system since the 
electron sheets around M is symmetric with respect to x-y exchange. These features, as well as the 
absence of Fermi surface nesting, suggest that the parent KFe2Se2 could be regarded as an electron 
doped 11 system with possible local moment magnetism. 

PACS numbers: 74.70.-b, 74.25.Ha, 74.25.Jb, 74.25.Kc 
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The discovery of iron-based compounds, typically rep- 
resented by LaFeAsO^. (1111-type), BaFejAsa^ (122- 
type) and FeSe'^'^(ll-type), has triggered enormous en- 
thusiasm in searching for the new high transition tem- 
perature superconductors without copper—^—. Ever since 
the discovery, density functional studies have been per- 
formed to explore the electronic structure and the pairing 
mechanism of the system. Calculations have been per- 
formed on LaFeAsOiSrJ^, and the ground state is found 
to be a coUinear anti-ferromagnetism (COL) state. The 
magnetic ordering was proposed to be the consequence 
of the Fermi surface nesting phenomen a^"' ^ ^1 , which is 
also present in the BaFe2As2 parent compound^-. The 
Fermi surface nesting is thus considered closely related 
with the superconducting (SC) phenomena since it is 
suppressed in SC phase. On the other hand, the J1-J2 
Heisenberg model based on a local moment picturei^ii^ii^ 
was also proposed to account for the magnetic struc- 
ture. It is also found that the band energy dispersion 
of these compounds should be calculated using the un- 
relaxed experimental structure in order to compare with 
the experiment o^^i^" , and that the ordered magnetic mo- 
ments on Fe atom are systematically overestimated in 
density functional calculations. 

Recently, another substance with similar chemical 
composition, KFe2Se2, has been produced^!. The mate- 
rial is reported to be iso-structural to BaFe2As2, and the 
superconductivity as high as 30 K is reported when it is 
intrinsically doped (Ki_a:Fe2Se2 with x — 0.2). Question 
thus arises that whether this material represents a new 
family or it is one of the discovered class. More specifi- 
cally, since KFe2Se2 is structurally close to BaFe2As2 but 
chemically close to FeSe, it is interesting to clarify which 
one is closer to its electronic structure. 

In this paper, we report our first-principles study of 
this compound. We demonstrate that the electronic 
structure of the parent KFe2Se2 could be regarded as 
an electron doped 11 system, instead of the structurally 
much closer BaFe2As2. All the calculations were per- 



formed with the Quantum ESPRESSO code^^, while 
an accurate set of PAW data^- were employed through- 
out the calculation. A 48 Ry energy cut-off ensures the 
calculations converge to 0.1 mRy, and all structures were 
optimized until forces on individual atoms were smaller 
than 0.1 niRy/bohr and external pressure less than 0.5 
kbar. For non-magnetic (NM) and checkerboard anti- 
ferromagnetic (CBD) states, a 8 x 8 x 8 Monkhorst- 
Pack k-grid^i was found to be sufficient; while for the 
coUinear anti-ferromagnetic (COL) state and bi-collinear 
anti-ferromagnetic (BIC) state, 6x6x8 and 16 x 8 x 4 
Monkhorst-Pack k-grid were needed to ensure the conver- 
gence to < 1 meV/Fe, respectively. The PBE flavor of 
general gradient approximation (GGA) to the exchange- 
correlation functionalS§. was applied throughout the cal- 
culations. 



We first examine several possible spin configurations 
for KFe2Se2 (TABLE |T|. The column expt indicates cal- 
culation with experimental structure and non-magnetic 
spin configuration, while the structures are fully op- 
timized (lattice constants as well as internal coordi- 
nates) for NM/CBD/COL/BIC columns. It is there- 
fore apparent that the coUinear phase, which is 57 
meV/Fe lower than the NM phase, is the ground state of 
KFe2Se2. This magnetic state ordering was also double- 
checked with full-potential linearized augmented plane 
wave (FLAPW) method using the elk code^. A body- 
centered tetragonal (bet) to base-centered orthorhombic 
(bco) phase transition is also present in the process, al- 
though the orthorhombicity e — 1 — b/a = 0.03% is al- 
most negligible. Although the phase transition is not 
yet observed in the experiments, the resistivity measure- 
ment shows an abrupt change around T — lOOK^^, which 
we propose to be due to the bco-bct phase transition. 
The magnetic moment on Fe atoms turns out to be 2.26 
/is/Fe for the coUinear phase, which is similar to those 
obtained in PAW calculations for BaFe2As2^. In order 
to estimate the magnetic coupling strength, we incorpo- 
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rate the J1-J2 Heisenberg model, defined by 

H = Ji ^ Si ■ Sj + J2 ^ Si ■ Sj 



<<ij>> 



Where, Si is the spin operator (of magnitude S) at the 
site i, < i,j > and << i^j >> denote the summation 
over the nearest neighbor and the next nearest neighbor 
sites, Ji and J2 are the nearest neighbor and the next 
nearest neighbor exchange interactions, respectively. By 
calculating the total energy per Fe atom for the CBD and 
COL states and assuming S = 1, we obtain J2 = 29 meV 
and Ji = 38 meV, respectively. 



TABLE I. Geometry, energetic and magnetic properties 
of KFe2Se2. Results in column expt were obtained us- 
ing the experimental structure and spin-unpolarized calcu- 
lations; while the NM/CBD/COL/BIC columns correspond 
to non- magnetic/checkerboard AFM/coUinear AFM/bi- 
collinear AFM configurations using the fully optimized struc- 
tures (lattice constants as well as internal coordinates), re- 
spectively. Ea is the total energy difference per iron atom 
referenced to the fully optimized NM structure, and mpc is 
the local magnetic moment on Fe. 





expt 


NM 


CBD 


COL 


BIC 


a (A) 


3.9136 


3.8791 


3.9058 


5.5930 


3.8271 


b(A) 


3.9136 


3.8791 


3.9058 


5.5916 


7.9530 


c(A) 


14.0367 


13.4476 


13.6849 


13.8525 


14.4783 


Ea (meV/Fe) 


272 





-18 


-57 


-39 


mpc (/ib) 


0.0 


0.0 


1.49 


2.26 


2.58 



Calculations with LaFeAsO, BaFe2As2, and FeSe 
systems suggest that the band dispersions and DOS 
of these systems should be studied without struc- 
tural optimization in order to compare with the 
experiments^^'^'^'"'^-ii^'22r— ; although their energetic 
properties as well as magnetism should be explored with 
structural relaxation. We followed this procedure, and 
the discussion in the rest of this paper will primarily focus 
on the calculations with unrelaxed (experimental) struc- 
ture unless we explicitly specify. Firstly, we present the 
density of states (DOS), as well as the projected density 
of states (PDOS) calculations (FIG. [J). The DOS and 
PDOS of KFe2Se2 resemble those of BaFe2As2 systems, 
and exhibit typical characteristics of the layered struc- 
ture. The contribution from Fe-3d and Se-4p orbitals 
dominates the states from Ep — 6.0 eV to Ep + 2.0 eV, 
while most of the K-4s contribution locates from Ep + 2.0 
eV to Ep + 6.0 eV. A closer examination of the PDOS 
data shows that over 90% of the states from Ep — 2.0 
eV to Ep are from the Fe-3d orbitals, and that the Fe- 
3d/Se-4p orbital hybridizes considerably from Ep — 6.0 
eV to Ep - 3.2 eV and from Ep to Ep + 2.0 eV. 

We further calculated the band structure for KFe2Se2, 
as shown in FIG. [2] Since the Se atom (4s^4p'') outer- 
most shell has 1 more electron than the As atom (4s^4p'^), 
the FeSe layer could be regarded as highly electron- 
doped. In fact, the band structure of KFe2Se2 shows that 
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FIG. 1. Total and projected density of states of KFe2Se2. The 
upper panel (solid line) is the total DOS; middle panel (dashed 
line) is the PDOS on Fe-3d orbitals; lower panel (dotted line) 
is the PDOS on Se-4p orbitals. We show only the energy 
range from Ep - 6.0 eV to Ep + 2.0 eV. 



the Fe-3dzx{y) and Fe-3d2.2_j,2 bands are fully occupied, 
whereas these bands were the origin of the hole pockets 
in the BaFe2As2 systems. At F, the 3dzx and 3dzy bands 
becomes degenerate due to the crystal symmetry. We de- 
fine the energy difference between these two bands and 



the Fermi level at F to be = Ep — ej 



and the differ- 



ence between these two bands and the 3da;2_y2 band at 



F to be As = £j 



(FIG. El). The latter is due 



to the slightly deformed tetrahedral crystal field by the 4 
Se atoms around the Fe atom. For the KFe2Se2 systems, 
et and As are 18 meV and 13 meV, respectively; while 
for the BaFe2As2 systems, they are -297 meV and 204 
meV, respectively. Similar to the BaFe2As2 band struc- 
ture, the bands close to Ep from T to Z are mostly flat, 
except for the one cross the Fermi level which is due to 
the hybridization of Fe-3d and Se-4p orbitals. It is wor- 
thy noting that the structural relaxation does not change 
the number of bands across the Fermi level, nor the or- 
bital character of these bands for KFe2Se2, in contrast 
to the cases for LaFeAsO and BaFe2As2. However, the 
structural optimization expands the band splittings A^ 
to 314 meV, and shifts the top of fully occupied d bands 
tt to 149 meV. 

The band structure we obtained is then fitted us- 
ing the maximally localized Wannier function (MLWF) 
method'^'^'"^^ (FIG. [2]) to obtain a model Hamiltonian for 
reconstructing the Fermi surfaces. To perform the fit- 
ting, we chose the 16 bands from Ep — 6.0 eV to Ep + 2.0 
eV, and 16 initial guess orbitals including the Fe-3d and 
Se-4p to ensure the fitting quality. Nevertheless, it is pos- 
sible to fit the band structure with slightly worse quality 
with the 10 Fe-3d orbitals only, in order to reduce the 
Hamiltonian size. Applying the symmetry, the number 
of orbitals could be further brought down to five. 

We show the parent KFe2Se2 Fermi surface in FIG. 
The Fermi surface of the parent KFe2Se2 consists 
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of two sheets around M-point and one sheet around F. 
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FIG. 2. Band structure of non-magnetic KFe2Se2 calculated 
with the experimental structure. The solid line is the den- 
sity functional theory (DFT) result and the dashed line is fit- 
ted with the maximally localized wannier function (MLWF) 
method. 



Despite of the similarities between the KFe2Se2 and the 
BaFe2As2 systems, two distinctions are apparent. First 
of all, the two sheets around M points in KFe2Se2 are 
much more symmetric than those in BaFe2 As2 . If we de- 
fine kp{6) to be the k- vector from the M point to the 
Fermi surface sheet, where 9 denotes the angle formed 
by the vector and M-M' (as shown in FIG. 3(c) I, we 

could further define t] — 1 — ,.,^_^/^Lo\ • For the two sheets 



fcF(-45°, 

around M in KFe2Se2, both yield r] < 1%, while for the 
BaFe2As2 system the inner sheet has rj — 48% and the 
outer sheet has ry = 41%. This feature suggests that the 
electronic structure of KFe2Se2 is in fact much closer to 
FeSe instead of BaFe2As2. Secondly, the sheets around T 
are completely different between KFe2Se2 and BaFe2As2 
or FeSe. Three sheets were observed for BaFe2As2 or 
FeSe system, which constitutes the three hole pockets 
for the system. For KFe2Se2 system, only one sheet ex- 
ists around (0,0, fcz) axis, which is highly 3-D and van- 
ishes around F. Thus, the Fermi surface nesting effect is 
absent in the parent KFe2Se2 compound. Nevertheless, 
one could achieve FeSe-like fermi surface using the rigid 
band model simply by shifting down the fermi level (FIG. 
|3(d)[ ) , or effectively by hole doping. From the band struc- 
ture calculation, the Fermi level has to be shifted down 
by 0.2 eV in order to recover the Fermi surface nesting 
effect. The DOS result indicates that shifting down Ep 
by 0.2 eV corresponds to 1.0 |e| hole doping effectively, or 
completely removing K from KFe2Se2. Due to the loss of 
Fermi surface nesting, the magnetism of KFe2Se2 is not 
simply due to the Fermi surface nesting effect. Instead, 
it is possible that the localized Fe-3d orbitals plays an 
essential role in the magnetism. Using the same model, 
we could also obtain the Fermi surface for Ki_2;Fe2Se2 
{x = 0.2), as shown in FIG. |3(b)[ which could be an 
analogue to an electron-doped 11 system. 

Finally, we study the {/-dependence of the magnetic 
coupling strength Ji and J2 using the GGA+U method. 
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FIG. 3. Fermi surface of Ki_2,Fe2Se2 reconstructed using the 
MLWFs. 
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3(c) for the sake 



3(a) and 3(c) are the plots for the parent 
compound {x = 0.0); panel 3(b) is the plot for x = 0.2; panel 
3(d) is obtained by shifting Ep by -0.2 eV. Panel 3(c) and 
3(d) are the 2-D plots of cross-section at fez — 0.0. Only the 
kp vector of the outer sheet is drawn in Fig. 
of visibility. 



to test if the COL configuration remains as the ground 
state if there is electron correlation beyond LDA. A se- 
ries of U from 1.0 eV to 6.0 eV were used to optimize 
the lattice constants as well as the internal coordinates 
for NM/CBD/COL configurations, and then Ji and J2 
under different U were fitted using mpe = LO/x^. Both 
Ji and J2 show linear dependence with respect to the on- 
site energy U, and the collinear state remains the ground 
state within a reasonable U range. 

In conclusion, we have studied the electronic struc- 
ture of KFe2Se2 using first-principles calculations. The 
ground state of KFe2Se2 turns out to be collinear anti- 
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ferromagnetic configuration with 2.26 ^,3 magnetic mo- 
ment on Fe atoms, similar to BaFe2As2. The Ji and 
J2 coupling strengths are calculated to be 0.038 eV and 
0.029 eV, respectively. Although the band structure is 
similar to heavily electron-doped BaFe2As2, the Fermi 
surface suggests that the system is much closer to an 
electron-doped FeSe system. The Fermi surface nest- 
ing effect is absent in the parent KFe2Se2 compound, 
thus the antiferromagnetism is possibly due to the lo- 
cal moments instead of the itinerant electrons. Since the 
KFe2Se2 is electron-doped, the superconductivity could 
be introduced with hole-doping or, effectively, K or Fe 
deficiencies. 
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Note added: After submitting this work to arXiv, 
we became aware of two recent papers^i^ where some 
related calculations have been also performed on the 
KFe2Se2 compound. 
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